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Abstract 

In dimension D < 2 the low temperature behavior of systems enjoy- 
ing a continuous symmetry is dominated by super-instantons: classical 
configurations of arbitrarily low energy. Perturbation theory in the 
background of a super-instanton produces thermodynamic answers for 
the invariant Green's functions that differ from the standard ones, but 
only in non-Abelian models and only starting at 0(1//?^). This effect 
modifies the /3-function of the 0{N) models and persists in the large N 
limit of the 0{N) models. 



Consider a Z lattice and let A be some given finite subset of it. The 
non-linear 0{N) a model with standard nearest neighbor interaction (s.n.n.i.) 
at inverse temperature (3 is defined by the following partition function Z: 

Z = [ Ylexp{f3S, ■ n('^(^' - 1)^^^) (1) 

Here S is an A^-vector, the integration is over the spins inside A with some 
boundary conditions (b.c). 

We shall be interested in the behavior of fixed (lattice) distance Green's 
functions, such as the expectation value of the energy density {S{0) ■ 5(1)), at 
large (3. As /3 — >■ cxd, this finite system will 'freeze' into the one or several con- 
figurations which minimize the energy (at given b.c). For jS large but finite, 
the spins will fluctuate mildly around one of these classical configurations and 
thus one can compute the Green's functions by a standard saddle point expan- 
sion. This perturbation theory (PT) approach produces the correct asymptotic 
expansion in more precisely, there exists an infinite set of finite numbers 
{dk} such that for /5 sufficiently large 

\{SiO)-Sil))-±^\<p, (2) 

Here {c^} are the expansion coefficients produced by PT. 

Of course both {c^} and {dk} depend upon the linear size L of A and 
the b.c. used. In particular, they could diverge as L — > oo. But it could also 
happen that even though the coefficients Ck have finite limits, for all k > ko the 
left hand side of eq.(2) is not o{(3~^) uniformly in L, i.e. if multiplied by (3^ it is 
not bounded by a function of (3 alone that vanishes in the limit (3 oo. This 
would indicate that the L — >■ oo limits of the Ck are not the coefficients of the 
correct asymptotic expansion in that limit. In other words, PT itself cannot 
be used to establish its correctness in producing the asymptotic expansion for 
(S'(O) ■ 5'(1)) at large (3. For that purpose one needs bounds uniform in L on 
the numbers {dk{L)}. This property of PT was established for the Abelian 
case 0(2) but the authors could not extend the proof to the non- Abelian 
case o\n) N >2. 

The general attitude in the community has been to ignore this issue. In 
fact as we tried to emphasize in the past Q there are some good reasons to 
expect trouble for D < 2. Firstly in those cases, the Mermin- Wagner theorem 
guarantees that the basic assumption underlying a saddle point approximation, 
the existence of a dominant configuration, is violated as L — >■ cx). Secondly in 
ID, where one can obtain the true answer using the transfer matrix formalism, 
one can see explicitely that for N > 2, the coefficient of 1//?^ in the expansion 
of (>S'(0) • >S'(1)) is finite yet incorrectly given by PT; that means that if one 
truncated the series at this order, the remainder would be 0(l//5^), rather 
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than o(l//3^), as required by a putative asymptotic expansion. 

In this letter we would like to point out that also in 2D, the two limits 
/5 — >■ cx) and L ^ oo do not commute for the non-Abelian case N > 2. The 
idea is the following: let A be a square of size L x L and on its edge freeze 
the spins along the x-axis. Since by the Mermin- Wagner theorem the infinite 
volume Gibbs state is 0{N) symmetric, for L — > oo, the expectation value of 
{S{0) ■S{1)) should be unchanged if the spin at the origin S{0) is also frozen in 
some given direction. So if one believes in PT, one should also expect that the 
PT coefficients approach values independent of whether S{0) has been frozen 
or not. As it will be seen, this is precisely what happens for 0(2) but not for 



First let us discuss the massless Gaussian, described by the partition func- 
tion 



It too possesses a continuous symmetry, which in 2D cannot be broken spon- 
taneously. Consequently, if we fix 99 = at the edge of the square and measure 
an invariant observable, such as -D(0, 1) = ((/^(O) — (/^(l))^, for L — cxd the result 
should become independent of whether V'(O) is fixed or not. In particular even 
if ip{0) = 0, D{0, 1) should converge to 1/2 for L — >• 00. 

This is clearly the case, as it can be seen from the data in Tab.l (line 1). 
These data were produced by a Monte Carlo study of this model using the over- 
relaxation procedure. Besides D{0,1) we measured certain Green's functions 
of the massless Gaussian with these b.c, which will be needed to compute 
the PT expansion of {S{0) ■ S{1)) in 0{N) - see below. We increased the 
number of measurements with L; as an example, for L = 80 we took 100, 000 
measurements, separated by 5 sweeps of the lattice. Similarly, the number 
of thermalization sweeps increased with L, reaching 10, 000 for L — 80. We 
estimated the error by a standard binning procedure. 

Let us return now to the case of interest of the 0{N) models. Firstly, 
although it is obvious, we verified that as we increase L, {S{0) ■ S{1)) does 
become independent of what we do with S{0). We studied both 0{2) and 
0(3). The spins at the edge of the square were fixed along the x-axis. The 
spin at the origin S{0) was either free ('Dirichlet b.c.') or also frozen along 
the X-axis ('super-instanton b.c.'). We performed a Monte Carlo study using 
a cluster algorithm. Some typical results are shown in Tab. 2, where we give 
also the value of {S{0) • S{1)) at the same /3 using periodic b.c. The results 
indicate that (>S'(0) • >S'(1)) does become independent of the b.c. as L — > 00, 
although the approach is slower for the super-instanton b.c. 

Next let us see what PT would predict for super-instanton b.c. Using the 
standard parametrization 



0(N) N>2 starting at 0{l/p^). 




(3) 
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5, = (vr„^l-7r2), 
then expanding in vr, one obtains 



(4) 



(5(0) -5(1)) 



where 



213 



+ (iV-l) 



G(l,l) 
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(5) 



G{^,J] 



{n{i) ■ 7r(j)), 

Y: [Gil,^r-Gil,JnG{^,^)-G{J,J)l 



<i,3> 



Y: [G{1,z) - G{l,j)]\G{l,j)[G{t,j) - G{j,j] 

<i,j> ^ 

Gil,t)[Git,j)-Git,t)] 



(6) 



The super-instanton b.c. on S* (vr = along the edge and at 0), imply that 
G{i,j) in eq.(6) is precisely the Green's function of the massless Gaussian 
studied above. From that study we already know that the term of 0(l//5) will 
converge to a value independent of the b.c. as L ^ oo. To study the term of 
0(1/ jS"^) we computed t and v numerically and report the values obtained in 
Tab.l. 

To interprete these data, let us recall that if instead of the super-instanton 
b.c. one would use periodic b.c, one would obtain M 



{S{0) ■ Sil)) = 1 



-D(l) 



N-1 



-D(lY + ... 



(7) 



f3 ^ 2/52 

Here D{1) = 2I) ~ x^) have neglected a term of 0(l//?2) coming from 

the elimination of the zero mode which vanishes very rapidly with increasing 
L. Since eq.(5) contains a term proportional to (A^ — 1)^, absent in eq.(7), 
the two would agree only if t converged to 1/8. In fact our data presented in 
Tab.l suggest that t converges to 0.035 or less. On the contrary, agreement of 
eqs.(5) and (7) for = 2 requires that v — t/2 converges to 1/8 and our data 
certainly seem consistent with that. 

The fact that v — t/2 goes to 1/8 can also be proven analytically. It follows 
after some algebraic manipulations that 
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v-tl1= G(0,z)-G(l,z)-G(0,j)+G(l,j) 

(8) 

It can be justified to take the limit L oo under the sum, and then it is 
seen easily by a summation by parts that it becomes 1/8. On the other hand, 
the same procedure is not justified for the sums defining v and t separately 
(because the sums are not convergent uniformly in L); one would obtain t = 
1/8 and v = 3/16, which would agree with Hasenfratz's formula (7), but clearly 
disagrees with our data. 

Therefore our results suggest that if in 2D one takes the limit L ^ oo 
termwise in the PT expansion, the result depends upon the b.c. used, but only 
for non-Abelian models and then only starting at 0(1/ (3"^). Since the correct 
asymptotic expansion has to be independent of b.c. by the Mermin- Wagner 
theorem, one must conclude that taking L ^ oo termwise is not a legitimate 
procedure. Of course it could happen that in 2D this procedure gives the 
correct infinite volume answer for certain b.c, just as it does in ID if one uses 
free b.c. For instance it has been argued that periodic b.c. must be giving 
the correct infinite volume answer because they agree with the 1/A^ expansion. 
Unfortunately, as we explained before 0, this argument is not correct because 
it involves again an a priori illegitimate interchange of two limits: /? — oo and 
N — > oo. So while it may happen that in 2D periodic b.c. do give the correct 
infinite volume PT coefficients, to our knowledge there is no support for that 
belief at the present time. 

We also computed the PT coefficients to the 2-point function at distances 
X larger than one lattice unit, and found that with super-instanton b.c. one 
obtains a different answer in the thermodynamic limit for N > 2. Since for 
large x {x > 4) this modification grows with x like it will give rise 

to a /5-function whose leading coefficient is different from the one found in the 
literature (obtained by using PT with periodic b.c). This fact casts doubt 
upon the alleged asymptotic freedom of the 0{N) N > 2 models. Details will 
be reported elsewhere. 

Over the years it has been argued that in theories which are perturbatively 
asymptotically free, PT should be applicable at short distances, where the sys- 
tem is well ordered. In fact starting at 0(1//?^), as it can be seen from eqs.(5) 
and (6), a PT computation involves sums over all distances. Consequently 
our finding that the trouble starts only from 0(1//?^) is quite understandable. 
The trouble does not arise in 0(2) because there, by parametrizing the spin 
as (cos(v9), sm{{p)), the Gibbs measure becomes a function of Vy? only. 

Before concluding let us explain why we called these b.c. super-instanton 
b.c. Let A be (the closest lattice domain to) a disk of radius R. Let us use 
the following spin parametrization 
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1 - TTfcOs{ipi)) 



(9) 



where the tTj are now (A^ — 2)-vectors. We choose as b.c. tt = along the edge 
and at 0, while ip{0) = but ip = c> along the edge. To do PT we need first 
the configuration of minimal energy with these b.c: it is some lattice version 
of the following continuum expression 



Please notice that this classical configuration has a vanishing energy as i? — >■ 
oo, hence the name super-instanton. 

More precisely, twenty years ago Kosterlitz and Thouless (KT) p argued 
that at large f3 the partition function is dominated by the classical configura- 
tions, which have since been called instantons. According to KT in 0(2) the 
dominant configurations are vortices and their energy is actually 0{log{R)) 
(hence they are bound), while in 0(3) they are 'hedgehogs' of energy 0(i?°), 
which create exponential decay no matter how large jS is. In fact, as we see 
now, the super- instantons are the dominant configurations in all 0{N) models, 
since their energy vanishes for R —y oo. They occur for all D < 2 and are the 
'enforcers' of the Mermin- Wagner theorem (a related idea underlies Pfister's 
[^] proof of that theorem): Let us for simplicity work in the continuum, with 
a playing the role of an UV cutoff, and use the fact that the configuration of 
minimal energy has a (f that only depends on the radius r. Then 



clogr/ logi? 







(10) 



E 



c 




(11) 



With the scalar product 




(12) 



we obtain 



c = (r 



l-D 



d_ 

' dr 



0) 



(13) 



and by Schwarz's inequality 




(14) 
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This gives a lower bound for E that goes to as i? oo, whereas for D > 2 
it remains bounded away from 0. Therefore one can have super-instantons 
(£; ^ yet c > Co > 0) only for D < 2. 

The arguments above suggest that at low temperature, the typical con- 
figuration in the 0{N) model in dimension D < 2 is a gas (or liquid) of 
super-instantons. The remarkable fact is that we had reached the same con- 
clusion two years ago, when we investigated these models by mapping them 
into a certain percolation problem ||^, ||; namely if we consider a certain patch 
of the sphere 5*^"^ of angular opening 0{l/^/J3) and ask what is its inverse 
image in a typical configuration, our percolation arguments indicated that the 
answer must be that it forms rings of arbitrarily large size (more precisely, 
neither the inverse image of the patch nor of its complement percolates). This 
super-instanton like structure was found to be responsible for the absence of 
exponential decay in all 0{N) models at (3 sufficiently large. Now we see that 
it may also be responsible for the failure of PT in non-Abelian models. Indeed 
imagine doing PT in the background of a super-instanton. Even though the 
new couplings induced by this external field vanish as l/(logi?)^, the potential 
IR divergences present in PT could produce nonvanishing effects as i? — > oo. 
In fact the effect could happen only for N > 2 {N = 2 involves only V0, 
which is IR finite) and only starting at 0(1//?^) since we must compensate 
two powers of 1/ logi? - this is in agreement with the findings reported above, 
which represented PT in the background of the trivial super-instanton. 

A. P. would like to acknowledge the hospitality of the IPN of the Universite 
de Paris Sud and of the Werner-Heisenberg-Institute. 
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Tab.l: The PT coefficients to order 1/(3'^ of{S(0) ■ S(l)) in the s.n.n.i. 0{N) 
models with super-instanton b.c. 



L 


10 


20 


10 


60 


80 


G(l,l) 


.3752(10) 


.4006(12) 


.4143(8) 


.4211(13) 


.4314(22) 


8?; - At 


.657(32) 


.709(40) 


.772(18) 


.792(18) 


.833(74) 


8t 


.345(46) 


.296(47) 


.305(18) 


.284(26) 


.317(77) 



Tab. 2a: {S{0)- S{1)) for the s.n.n.i. 0(2) model with periodic (per.), Dirichlet 
(Dir.) and super-instanton (s.i.) b.c. (Monte Carlo data taken at j3 = 1.1; the 
data point for L = 256 is taken from R.Gupta and C.Baillie, Phys.Rev. B45 
(1992) 2883.) 



L 


20 


40 


80 


256 


per. 


.7177(6) 


.7155(3) 


.7146(2) 


.7138(1) 


Dir. 


.717(2) 


.712(2) 






s.i. 


.774(3) 


.768(3) 


.755(4) 





Tab. 2b: (S'(O) ■iS'(l)) for the s.n.n.i. 0(3) model with periodic (per.), Dirichlet 
(Dir.) and super-instanton (s.i.) b.c. (Monte Carlo data taken at P = 2.0; 
the data point for L = 1024 is taken from J.Apostolakis et al, Phys.Rev. D43 
(1991) 2687.) 



L 


20 


40 


80 


1024 


per. 


.7272(5) 


.7257(3) 


.7253(2) 


.72511(3) 


Dir. 


.728(2) 


.726(2) 






s.i. 


.776(3) 


.766(3) 


.753(4) 
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